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Abstract 



We present a complete description of a new physical model with goldstino supermultiplet for Com- 
pHEP software package. The model contains MSSM content supplemented by three new fields: goldstino 
- two component spinor, scalar sgoldstino and pseudoscalar sgoldstino. To the leading order in 1/F, 
where F is the parameter of supersymmetry breaking, and to zero order in MSSM gauge and Yukawa 
coupling constants, the interactions between the component fields of goldstino supermultiplet and MSSM 
fields are derived. They correspond to the most attractive for collider studies processes where only one 
of these new particles appears in a final state. In this case light gravitino behaves exactly as goldstino. 
For sgoldstinos, as they are R-even, we include in the package only sgoldstino couplings to goldstino and 
sgoldstino couplings to SM fields as the most interesting phenomenologically. All new coupling constants 
between the components of goldstino supermultiplet and MSSM fields are completely determined by 
' the ratios of soft supersymmetry breaking parameters and F. Thus, in addition to the MSSM coupling 

constants there are three new parameters in the package: the scale of supersymmetry breaking, and 
^ ^ masses of scalar and pseudoscalar sgoldstinos. The CompHEP notations for the new fields and parame- 

D , ters are presented. The package was successfully tested and all relevant results known in literature have 

• been reproduced with good accuracy. It may be used in order to calculate any tree-level process with 

one gravitino or sgoldstino on mass shell. Universality of the Lagrangian allows to apply the package in 
study of phenomenology of any supersymmetric extension of the Standard Model. 

X 

H , 

5^ , 1 Introduction 

CompHEP package is adjusted to calculate automatically particle decay rates and cross sections at the 
tree level.|^ It is aimed to improve the accuracy, to cut down the efforts and to shorten the time usually 
required in study of the collision processes at high energy. 

First versions of CompHEP contained physical models based only on the set of particles and interactions 
forming the Standard Model (SM). In CompHEP any user has the opportunity to change vertices and add 
new particles into the SM content. Thereby some modified models (based on the SM) have been created 
and used in calculations (see e.g. Refs. |^ for leptoquark models and Refs. ||] for SM with anomalous 
couplings). The growth of activity in a study of phenomenology of various extensions of the SM inspired the 
authors of the package to include the Minimal Supersymmetric extension of the Standard Model (MSSM) 
into CompHEP^ Below we will refer to this model as CompHEP/SUSY one. 

In the framework of supersymmetry there are several types of models which give significantly different 
predictions for the values of low-energy MSSM parameters and for low-energy spectrum, that might be 
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tested at accelerators. The incorporation of different classes of models into CompHEP gives an opportunity 
to study quantitatively high energy phenomenology of these models getting results with high accuracy and 
with minimum boring routines. The first step to this direction had been done, when the possibility to 
work with the set of parameters evaluated in the frameworks of mSUGRA and GMSB became available in 
CompHEP/SUSY. The incorporation of other extensions, in particular, the MSSM with R-parity violation, 
are in progress. 

This paper is devoted to the description of the supersymmetric extensions of the SM with very light 
gravitino and light sgoldstinos. These new particles as well as their interactions with MSSM fields have been 
included in the new version of CompHEP/SUSY model.^ The paper is organized as follows. In section 
2 we describe the supersymmetric model with sgoldstino supermultiplet. Namely, section 2.1 contains the 
description of goldstino supermultiplet. In sections 2.2 and 2.3 we present the effective Lagrangians for 
very light gravitino and light sgoldstinos, respectively. In section 2.4 we derive the free Lagrangian for 
light gravitino and sgoldstino fields, as well as the effective interactions responsible for sgoldstino decays 
into gravitinos. A sketch of phenomenology of the models with very light gravitino and models with light 
sgoldstinos is presented in Section 3. It is aimed to outline the experimentally allowed region of parameter 
space of the model. Section 4 explains the CompHEP notations for the new particles and parameters 
additional to ones used in CompHEP/SUSY model. Section 5 is devoted to concluding remarks, the tests 
which have been done to check the incorporated models, and possible ways to improve the package. 



2 Description of the model 
2.1 Goldstino supermultiplet 

In any supersymmetric extension of the Standard Model (SM) of particle physics spontaneous supersymmetry 
breaking occurs due to non-zero vacuum expectation value of an auxiliary component of some chiral or vector 
superfield. As the most simple case, let us consider a model where the only chiral superfield 

S = s + V20'tjj + 9^F, (1) 

obtains non-zero vacuum expectation value F for the auxiliary component, 

(^.> ^F. (2) 

Then ip is two-component Goldstone fermion, goldstino, and its superpartners 

S^±{s + sn, P^-^{s~s*), (3) 

are scalar and pseudoscalar sgoldstinos, respectively. 

In the framework of supergravity dfj^tp becomes a longitudinal component of gravitino, that is the super- 
Higgs effect. As a result, gravitino gains mass m^^2 which is completely determined, in a realistic model 
with zero cosmological constant, by the parameter of supersymmetry breaking F: 

VS^ F 

^3/2 = ^^7 ■ (4) 

Vo Mpiank 

Sgoldstinos remain massless at tree level and become massive due to corrections from high order terms in 
the Kahler potential. Provided these terms are sufficiently suppressed sgoldstinos are light and may appear 
in particle collisions at high energy colliders. Such pattern emerges in a number of non-minimal supergravity 
models Q and also in gauge mediation models if supersymmetry is broken via non-trivial superpotential 
(see, Ref. |^ and references therein). Here we consider sgoldstino masses as free parameters. 

Gravitino and sgoldstinos interact with MSSM fields and corresponding coupling constants are inversely 
proportional to the parameter of supersymmetry breaking F . This implies that gravitino has to be very light 
(see Eq. (H)), otherwise gravitino as well as sgoldstinos is effectively decoupled from the MSSM fields at the 
energy scale which is expected to be available for study at the nearest future colliders. 
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2.2 Gravitino effective Lagrangian 

The effective Lagrangian for the gravitino is obtained from the N=l supergravity theory and may 
be used to calculate scattering processes involving all helicity components of the massive gravitino. Mean- 
while the energy scale attainable at the present and the nearest future generation of accelerators favors the 
longitudinal component of gravitino as the most promising to be searched for in collision experiments. The 
matter is that while the effective coupling constants for ±| helicity components of gravitino are suppressed 
by gravity scale, the longitudinal component of gravitino, in accordance with the Equivalence Theorem [0, 
effectively behaves as a two-component fermion, goldstino %p, 

G^ ^ i df^il; , (5) 

v3m3/2 

so gets enhancement if the energy E of the process is significantly larger than gravitino mass. Therefore, 
goldstino couplings are suppressed by the ratio E/F (see Eq. (^)), which may be significantly lower than 
gravity scale. The effective low energy Lagrangian is non-renormalizable, and on general grounds the sensi- 
tivity of collider experiments to the scale of supersymmetry breaking VF is expected to be not better than 
about 10 TeV. Thus, in the models under discussion gravitino is very light, 7713/2 ^ a few ■ 10~^ eV (lighter 
than a few tens millielectronvolts), and only longitudinal component, goldstino, is relevant. Hereby gravitino 
is considered as massless {Mpiank —>■ 00, F = fixed) in this approach. 

Further simplification is related to the fact that cross sections of the processes with few goldstinos in 
a final state are suppressed by additional powers of F in comparison with single goldstino production. 
Consequently, searches for the latter process are expected to be more sensitive to the scale of supersymmetry 
breaking, so we retain only one-goldstino couplings in the effective Lagrangian. In this case the interaction 
between goldstino and other fields is plainly given by Goldberger-Treiman formula 



with JgijgY being a supercurrcnt. For MSSM this interaction reads 

V all gauge 

fields 



F 



all matt 
Higgs field 



where a and k run over component fields of vector (A", v" ^) and chiral {(fik, <Pk) MSSM supermultiplets, 
respectively (see Eq. (||) for details); a'^ — (!,(?), (t^ = (1,— (t). Throughout this paper we adopt two- 
component fermions. 

Whereas we are interested in processes with light gravitino at an external line of Feynman graphs (i.e., 
tree-level production of gravitino), we simplify Eq. (|6|) by making use of the goldstino equation of motion 
(equation of motion of a free massless two-component spinor field). Finally we obtain 



■ d Higgs fields 



Such interaction terms exist in any supersymmetric extension of the SM, so this Lagrangian may be regarded 
as the universal tool to study the sensitivity of various processes to the scale of supersymmetry breaking. 

It is worth noting, that as the total amplitudes with helicity ±i gravitino in a final state have to 
vanish unless supersymmetry is broken, these amplitudes should be proportional to some supersymmetry 
breaking parameters. Therefore one can use the equations of motion of MSSM fields to rewrite the effective 
Lagrangian (^ in the form, where all vertices are proportional to soft breaking terms. Although this form is 
very eligible for simple estimates, since the leading order dependence of the amplitudes on the soft parameters 
is explicit, CompHEP operates with the Lagrangian (M) because of some technical reasons. 
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2.3 Sgoldstino effective Lagrangian 



We obtain the low-energy effective Lagrangian for sgoldstino by making use of spurion method g . It consists 
in exploiting the fact that, by definition, sgoldstinos are scalar components of the very supermultiplet S whose 
auxiliary component gains non-zero vacuum expectation value due to supersymmetry breaking. Then one 
can consider a simple supersymmetric model with non-renormalizable interactions between S and MSSM 
superfields, which provide MSSM soft terms if non-zero (Fg) is taken into account. Consequently, the 
corresponding coupling constants are fixed by the ratios of the soft terms and F. 

We begin with the Lagrangian of the MSSM, which consists of two parts: the first part is a minimal super- 
symmetric generalization of the Standard Model (SM), while the second one represents the supersymmetry 
breaking, 

^MSSM = ^SUSY + ^breaking , 

where 



^SUSY = ^gaii 



hie 



■' superpotential 



and 



r 



gauge 



d^e TtW^W" + h.c. 



hie 



E 



all matter 
id H iggs fields 



c 



superpotential 



where L, E and Q, U, D are lepton and quark superfields, respectively, Hd and Hu are electroweak Higgs 
superfields; i,j = 1, 2 are SU(2) and a,b — 1, 2, 3 are generation indices, color indices are omitted and two- 
component antisymmetric tensor is defined as ei2 = — 1; yab are matrices of Yukawa constants and is a 
parameter of higgsino mixing. Hereafter we assume that R-parity is conserved. The expansion of superfields 
in series of grassmanian variables 9,9 introduces the component fields: 

= 0fe + V29(j)k + e^Fk , 

V = T'^Va , Tr T-^T^ = S'^'' , K = -9a^'9va ^ + i9W\a + ^WD^ , 
W = T'^W: , W: = -iK + OD<i ~ '-a^a'^9F^^, + 9^<j^^dX , 

K = . - + fabe Vt^v^ , , D^\l ^d^\l+ fabe . («) 

where for vector superfields we chose the Wess-Zumino gauge. 
Generally, supersymmetry breaking part has the form 

-Areafc^ng =^m2|0fe|2+ E TrA"A" + /l.C.^ 

- 6,, (Bh^ohlj + A^JiilWj, + A^MWj, -f A'i.it.Ulhlj + h.c) , (9) 

where k (a) runs over all scalar (pk (gaugino Aq,) fields. 

Since supersymmetry is broken spontaneously, Eq. (j^) implies the following effective interaction between 
MSSM superfields and goldstino supermultiplet S: 



C-S-MSSM — Cs-Kdhler + ^S- 



+ Cs-superpotenUal , 



where 



Cs-Kdhler = - ld?9 CP9 S ■ ^ 



II matter 
Higgs fields 



.2 
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h.c. 



all gauge 
fields 

f / B A^ \ 

CS-superpotenUal = J<feS- HhEfj + LiE^.H}, + QiD^.H}, + QlU^H^u) + ^-C 

These terms emerge if the fields from the hidden sector, where supersymmetry breaking occurs, are integrated 
out. The only remnant of the hidden sector is goldstino supermultiplet, which is supposed to be light. 

Integrating over 0, 9 and taking into account Eqs. (0), (H) we obtain soft supersymmetry breaking 
terms as well as interaction between the components of goldstino supermultiplet and components of 
MSSM superfields: 



•'S—K'dhler 



all matter 
id Higgs fields 



E {^■s-^lF, + h.c.) , (10) 



fillets' 



-S-superpotential — —^ij ( ' XdXu ^ ■ i^^jjFfj^ + ^Ho^U 



1 



(9>h • shlj + sql ■ utxh + sul ■ qlxh ~ ^F^^AK " sqlFu^hlj - sq'^ulFjj^ 



(12) 

Here we presented only leading order in ^ terms; the convention for the Levi-Civita tensor is e*'^^^ = — 1. 

Eliminating auxiliary fields, one gains the low-energy effective Lagrangian for the interaction between 
the components of goldstino supermultiplets and MSSM fields. A few remarks are in order. First, being 
aimed at getting the most universal Lagrangian applicable to any supersymmetric extension of the SM, we 
have to consider only interaction terms obtained to the leading order in and to zero order in MSSM 
coupling constants. Second, from the phenomenological point of view the most promising interaction terms 
are ones which include the minimum new (additional to SM) fields. In model with R-parity goldstino is 
R-odd, whereas sgoldstinos are R-even. Consequently, the simplest vertex with goldstino has to contain 
at least one new field - superpartner, while among sgoldstino vertices there is a set containing only SM 
fields, and only these vertices are of interest. Likewise we reject the interactions with F and D terms, since, 
being resolved, these auxiliary fields are proportional to SM gauge or Yukawa couplings. Thus, only minimal 
sgoldstino-partner-partner vertices have been included in the presented model: 



E 

all gauge 
fields 



4V2F " 



/ aL aD aU 

- [-^^i^t ■ Sh], + ^qidt ■ Shh + ^9>6 • Shi ] ] + h.c. , (13) 

(M 
E ^p-Pa\.e'^'"'F:,^ 

fields 
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/ AL aD aU \\ 

- ■ Phh + ^^lidt ■ Ph^D + *;^9>g • Ph'uj j + h.c. (14) 

Here we used the replacement (||). 

Being included into CompHEP, this model should be regarded as an additional option allowed in the 
framework of CompHEP/SUSY model. Currently accessible versions of this package operate only with real 
parameters and coupling constants. Likewise the trilinear soft couplings are assumed to be proportional to 
the corresponding Yukawa couplings, Aab = AabUab- The sgoldstino Lagrangian transformed in accordance 
with these rules reads 

M 



ad aU 

f^yf, • S{e.,qidlh}, + h.c.) - -^y^, ■ S{e,,<f,ulK + h-^) > (15) 



hi 



fiLids 

- ^^Vl ■ P^^^^aAlh^D - h.c.) - ^^Vab ■ P{e^,q>thlJ - h.c.) . (16) 

One can check that resulting interaction terms between sgoldstinos and SM gauge bosons coincide with 
corresponding couplings presented in Ref. except of the terms proportional to the parameter /i^ introduced 
there. The resulting interaction terms between sgoldstinos and SM fermions coincide with corresponding 
couplings presented in Ref. pOd up to overall sign. 



2.4 Gravitino— sgoldstino effective Lagrangian 

To obtain the effective Lagrangian for gravitino-sgoldstino sector we again apply the spurion method 
It is suitable for our purposes, since the only process we want to describe is the decay of sgoldstinos into 
gravitinos, i.e., the process with all involved particles being on mass shell. These decays are kinematically 
allowed in the region of parameter space viable for study at colliders. The relevant terms follow from the 
effective Lagrangian 

Ls = / dO^S^ (S^S - (5t5)^) - 4^ (/ dO^S^ + /..c.) , (17) 

where we introduced sgoldstino masses. Recalling the definitions (|l]), (^), (^ we obtain the Lagrangian 
with kinetic and mass terms for goldstino (longitudinal component of gravitino) and sgoldstinos and the 
interactions respectable for sgoldstino decays into gravitinos 

2 2 
' s{^i; + i;^)-t-P^P{i;^-^^) . (18) 



2V2F ^ ' 2V2F 

The results for sgoldstino-gravitino coupling constants coincide with corresponding terms presented in Ref. 
up to overall sign. 



3 Phenomenology of the model 

Till now, there are no experimental evidence for gravitino or sgoldstinos. The study of their phenomenology 
results in acquisition of bounds on their coupling constants. 

Let us briefly review the current status of experimental limits on F. 
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If we ignore sgoldstino (or, equivalently, set ms^p — > oo), then in models with light gravitino the current 
direct bound on the supersymmetry breaking scale is obtained from Tevatron, \/F >217 GeV fO]. The 



upgraded Tevatron may be able to cover the range of ^/F up to ^/F ~290 GeV |12 , while LHC will be 
capable to reach ^/F ~1.6 TeV As we will see, in models with light sgoldstinos, collider experiments 
become more sensitive to the scale of supersymmetry breaking. 

Then, in the model under discussion gravitino is extremely light: 111,3/2 ^ 10~^ eV and may affect 
standard Big Bang Nucleosynthesis as a neutrino-like specie. To get a sufficient dilution the gravitino 
decoupling temperature must exceed T ^ 200 MeV, that yields the limit on the scale of supersymmetry 
breaking, \/F > 600 GeV at mg ~ 100 GeV Similar bounds follow from the cooling rate of white 

dwarfs p5| and Z-boson width [ p^ in the models where one neutralino is lighter than Z-boson. It is worth 
noting that all these bounds are ensued from the study of the processes with two gravitinos, so they are 
model dependent]^ and may be invalid in other supersymmetric models. 

Let us turn to the models with light sgoldstinos. Note, that all sgoldstino coupling constants introduced 
in previous sections are completely determined by MSSM soft terms and the parameter of supersymmetry 
breaking F, but sgoldstino masses (mg, mp) remain free. Depending on the values of these unknown param- 
eters, sgoldstinos may show up in different experiments. Phenomenologically interesting models form four 
classes: 

1) Sgoldstino masses are of order of electroweak scale, y/F ~ 1 TeV — sgoldstino may be produced in 
collisions of high energy particles at colliders ^, |l^ . 

2) Sgoldstino masses ms, mp ^ 1 MeV-^1 GeV, VF ^ 1 TeV — sgoldstinos may emerge in products of 
rare decays of mesons 0, such as T — > S{P)j, J/ip — > S{P)j. 

3) Sgoldstinos in models with flavor violation in the sector of soft trilinear couplings, Aab 7^ ASab — 
sgoldstino interactions lead to flavor violating processes. Namely, sgoldstinos may contribute to FCNC (mass 
difference and/or CP-violation in the systems of neutral mesons) |l^, then, if kinematically allowed, they 
appear in the product of rare decays, such as t cS{P) ^ eS{P), K ttS [0, etc. 

4) Sgoldstinos are lighter than 1 MeV — these models may be tested in low energy experiments 0, such 
as reactor experiments, conversion in magnetic field, etc. Also sgoldstinos may play very important role in 
astrophysics and cosmology ^ |T^: they may change the predictions of Big Bang Nucleosynthesis, 
distort CMB spectrum, affect SN explosions and cooling rate of stars, etc. 

There are flavor-conserving and flavor-violating interactions of sgoldstino fields. 

As concerns flavor-conserving interactions, the strongest bounds arise from astrophysics and cosmology, 
that is Vf> 10^ GeV g |l§, or TO3/2 > 600 eV, for models with mgi^p) < 10 keV and MSSM soft flavor- 
conserving terms being of the order of electroweak scale. For the intermediate sgoldstino masses (up to a 
few MeV) constraints from the study of SN explosions and reactor experiments lead to VF > 300 TeV [Q . 
For heavier sgoldstinos, low energy processes (such as rare decays of mesons) provide limits at the level of 



F > 500 GeV [10 



The collider experiments exhibit the same level of sensitivity to light sgoldstinos as rare meson decays. 
Indeed the studies ^ |2^, ^ of the light sgoldstino {jns^p ^ a few MeV) phenomenology based on 
the effective low-energy Lagrangian derived from N=l linear supergravity yield the bounds: ^/F > 500 GeV 
(combined bound on Z ^ Sff^Pff [|6|; combined bound on e+er jS,-fP @) at M^oft ^100 GeV, 



VF > 1 TeV pTl (combined bound on pp gS, gP) at M3 ~ 500 GeV. Searches for heavier sgoldstino 
at colliders, though exploiting another technique, results in similar bounds on the scale of supersymmetry 
breaking. Most powerful among the operating machines, LEP and Tevatron, give a constraint of the order of 
1 TeV on supersymmetry breaking scale in models with light sgoldstinos. Indeed, the analysis carried out by 
DELPHI CoUaboration H yields the hmit y/F > 500-^200 GeV at sgoldstino masses ms^p = 10-=- 150 GeV 
and Msoft ~ 100 GeV. The constraint depends on the MSSM soft breaking parameters. In particular, it 
is stronger by about a hundred GeV in the model with degenerate gauginos. At Tevatron, a few events in 
pp — > S^{Z) channel, and about 10'* events in — > 5' channel would be produced at ^/F — 1 TeV and 
Msoft ~ 100 GeV for integrated luminosity C ~ 100 pb^^ and sgoldstino mass of order of 100 GeV |p7| . 
This gives rise to a possibility to detect sgoldstino, if it decays inside the detector into photons and \fF is 
not larger than 1.5 2 TeV. 



*A11 these bounds were obtained in the framework of hnear N=l supergravity. 
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The flavor-asymmetric processes are generally very sensitive to light sgoldstino. Namely, with flavor- 
changing ofF-diagonal entries in squark (slepton) squared mass matrix close to the current bounds, direct 
measurements of decays of mesons (leptons) provide very strong bounds, up to > 900(15000) TeV [ p^ 
(vahd at ms < 5(0.34) GeV), which is much higher than bounds expected from future colliders. Sgoldstinos 
may appear in top-quark decay, that could help to probe the scale of supersymmetry breaking at LHC up 
to about 10 TeV The possible contributions of sgoldstinos with ms,p = 100 GeV to mass differences 
in the systems of neutral mesons yield the bounds \/F > 10 TeV If off-diagonal entries are small, the 
limits on \/F become weaker: they scale as square root of the corresponding off-diagonal elements. 

Among the laboratory processes, the most sensitive to very light sgoldstinos are propagation of laser 
beam in magnetic field and reactor experiments. For heavier sgoldstinos the most sensitive processes to 
flavor-conserving sgoldstino couplings are T decays and associated sgoldstino production in collisions of 
high energy particles; charged kaon decays are most sensitive to flavor-violating sgoldstino couplings. For 
sgoldstino masses of the order of electroweak scale the most sensitive process to flavor-conserving sgoldstino 
couplings is sgoldstino resonant production at hadron colliders, while the most sensitive to flavor-violating 
couplings is measurement of mass differences in the systems of neutral mesons. 

4 CompHEP notations 

In the CompHEP package gravitino and sgoldstinos {ip, S, P) are denoted as ~g, sS and sP, respectively. 
Sgoldstino masses ms and mp are referred as MsS and MsP; gravitino is massless. The scale of supersymmetry 
breaking in the entire theory, Vf, corresponds to sqF. The Feynman rules for the Lagrangians (0), (^5|), 
( pj| ) , (|lj) have been obtained by making use of LanHEP package |^ . They have been incorporated into 
CompHEP, as well as aforementioned new parameters with default values MsS = MsP = 200 GeV, sqF = 
3 TeV.|^ Subsequent checks gave results consistent with ones known in literature. 

5 Concluding remarks 

We described the model with very light gravitino and light sgoldstinos, which has been included into Com- 
pHEP/SUSY package. This package may be used in order to calculate a rate of any tree- level process with 
one gravitino or sgoldstino being on mass shell. This includes: 

• decay of a superpartner into partner and gravitino 

• sgoldstino decays 

• gravitino inelastic scattering off MSSM particles into MSSM particles 

• sgoldstino inelastic scattering off SM particles into MSSM particles 

• gravitino (sgoldstino) production in collisions of MSSM particles. 

The effective Lagrangian incorporated into CompHEP/SUSY is universal, so the package may be used to 
study the phenomenology in any supersymmetric extension of the Standard Model. The only interaction 
terms to the leading order in 1/F and to zero order in MSSM gauge and Yukawa coupling constants, relevant 
for gravitino and sgoldstino phenomenology, are very terms included into this package. 

Note, that direct independent measurement of MSSM soft supersymmetry breaking terms and gravitino 
or/and sgoldstino couplings gives the unique possibility to estimate the scale of supersymmetry breaking 

Vf. 

As a flnal check for gravitino couplings we compared the analytical results (from CompHEP output) for all 
decay rates of the kind superpartner— ^partner + gravitino with corresponding formulae presented in Ref. 
and found no difference. To make sure that sgoldstino vertices are correct we compared the analytical results 
(from CompHEP output) for all sgoldstino decay rates with corresponding formulae presented in Refs. ^ 
The results coincide up to terms proportional to additional parameter /i^ introduced in Ref. To check 

^ These values may be changed by user that is a standard CompHEP menu option. 
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the numerical calculations carrying out by the package, we have reproduced with good accuracy all results 
for sgoldstino production in various channels presented in Refs. [T^ . 

There are several remarks which may be helpful for the user who are going to work with this package. 

It is worth recalling that we discussed the model with R-parity being conserved. Gravitino is R-odd, 
consequently, it might appear in collision of SM particles only accompanied by another R-odd particle. 
Moreover, gravitino is stable the lightest supersymmetric particle — LSP. Sgoldstinos are R-even, therefore, 
their production may be not accompanied by production of superpartners. Sgoldstinos may decay into SM 
particles, and since two-photon channel is always open, they are unstable. 

The second remark concerns the treatment of sgoldstino coupling constants to be proportional to the 
MSSM soft breaking terms. In fact, the Lagrangians (p^), ( p6| ) describe sgoldstino interactions at the 
superpartner scale. Sgoldstino coupling constants at other energies may be obtained by making use of 
renormalization group evolution. As an example, for the gluonic operator one has 

G^.(M3)=G^.(M)^^-^fi^ (19) 

where P{as) is QCD /3-function. These corrections (renormalization group enhancement) are neglected in 
the package, but may be taken into account plainly by making use of appropriate rescaling. 

Then, in CompHEP all Yukawa coupling constants are determined as the ratios of fermion pole masses 
to the corresponding Higgs vacuum expectation value. In order to estimate sgoldstino branching ratios into 
fermions one should take into account the difference between running Yukawa constants and Yukawas used 
in CompHEP. 

The forth remark is that sgoldstino couplings to superpartners are not included into this package. If open 
(allowed kinematically) , new sgoldstino decay channels distort the pattern of sgoldstino branching ratios 
calculated in CompHEP. In a given model the level of distortion depends on the set of soft supersymmetry 
breaking masses and sgoldstino masses and effective couplings. 

Finally, we discuss the possible directions to improve this package. The first interesting possibility is 
to take into account CP- and flavor-violating sgoldstino couplings. This option will be allowed, if complex 
coupling constants and flavor violating couplings become available in CompHEP. Then, sgoldstino couplings 
to superpartners may be also included. In the spurion approach corresponding coupling constants are not 
completely determined by MSSM soft breaking terms, but may depend on new parameters originated from 
the hidden sector. Similar situation takes place with sgoldstino interaction terms proportional to MSSM 



gauge or Yukawa couplings. For instance, the account of F and D auxiliary fields (see Eqs. (|10|), (|ll|), (|12[)) 
results, in particular, in sgoldstino-higgs mixing with coefficients depending also on new parameters from 
hidden sector. The model independence of these coefficients will be clarified elsewhere. 
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